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10 We investigate the chemical interaction between Saturn's corotating plasma 

11 and Enceladus' volcanic plumes. We evolve a parcel of ambient plasma as 

12 it passes through a prescribed H2O plume using a physical chemistry model 

13 adapted for water-group reactions. The flow fleld is assumed to be that of 

14 a plasma around an electrically-conducting obstacle centered on Enceladus 

15 and aligned with Saturn's magnetic field, consistent with Cassini magnetome- 

16 ter data. We explore the effects on the physical chemistry due to: (1) a small 

17 population of hot electrons; (2) a plasma flow decelerated in response to the 

18 pickup of fresh ions; (3) the source rate of neutral H2O. The model conflrms 

19 that charge exchange dominates the local chemistry and that H3O"'" dom- 

20 inates the water-group composition downstream of the Enceladus plumes. 

21 We also find that the amount of fresh pickup ions depends heavily on both 

22 the neutral source strength and on the presence of a persistent population 

23 of hot electrons. 
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1. Introduction 

24 Early Cassini encounters with Enceladus revealed surprising evidence of a significant 

25 source of water (with trace percentages of other neutrals, including CO2) from geysers 

26 located at the moon's southern pole [Hansen et al., 2006; Porco et al., 2006; Spencer 

27 et al., 2006; Waite et al., 2006]. The H2O cloud reacts with Saturn's corotating plasma 

28 torus, loading Saturn's magnetosphere with fresh ions. The pickup rate M quantifies 

29 the amount of fresh ions added to the magnetosphere from charge exchange and im- 

30 pact/photoionization. Respective contributions to pickup from charge exchange and im- 

31 pact/photoionization are fundamentally different in that charge exchange does not con- 

32 tribute to ion production because one ion replaces another. Both processes however 

33 introduce slow-moving ions which must subsequently be accelerated by Saturn's magne- 

34 tosphere. 

35 Early identification of the interaction between the local water source and Saturn's coro- 

36 tating plasma was made by Dougherty et al. [2006]. Based on the Cassini Plasma Spec- 

37 trometer (CAPS, Young et al. [2004]) analysis by Tokar et al. [2006], Pontius and Hill 

38 [2006] modeled the interaction and derived a pickup rate of M « 100 kg s~^. Khurana 

39 et al. [2007] and Saur et al. [2008] (hereafter K07 and SOS) discovered a range in M 

40 (0.2-3 kg s~^) from Cassini magnetometer data in the three earliest Cassini Enceladus 

41 flybys (EO, 17 February 2005; El, 09 March, 2005; E2, 14 July 2005). Constrained by Ion 

42 Neutral Mass Spectrometer (INMS, Waite et al. [2006]) and Ultraviolet Imaging Spectro- 

43 graph (UVIS, Hansen et al. [2006]) observations. Burger et al. [2007] estimated a pickup 

44 rate ol M ^ 2-3 kg s~^ from a neutral cloud model. The large discrepancy between 
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45 the pickup rates derived from CAPS and magnetometer data is due not only to the fact 

46 that the region considered by Pontius and Hill [2006] is much larger than that consid- 

47 ered by K07 and SOS, but also because the Pontius and Hill [2006] result depends on the 

48 poorly-constrained value of Saturn's Pederson conductivity. 

49 In this paper, we use a physical chemistry model to investigate the chemical interaction 

50 between the corotating plasma and the Enceladus plumes. Charge exchange dominates 

51 the local chemistry and leads to an HsO^-dominated plasma downstream of Enceladus. 

52 We find that pickup increases when hot electrons are present — more so with a high neutral 

53 source rate. 

2. Model 

54 We use a physical chemistry model developed to study the Enceladus torus [Delamere 

55 and Bagenal, 2003; Delamere et al., 2007; Fleshman et ai, 2009] to investigate the com- 

56 position of plasma traveling along prescribed flow lines. The model evaluates mass and 

57 energy rate equations for water-group ions (W"*" = O"*"-!- OH"*"-!- H20"'"-|- HaO"*"), pro- 

58 tons, and thermal electrons in a parcel of plasma transiting the simulation. Neutrals 

59 are assumed to be cold, and in this study neutral abundances are fixed. The full set of 

60 reactions includes charge exchange, photoionization, ionization by electron impact, radia- 

61 five excitation, recombination, and molecular dissociation by both electron impact and 

62 recombination. All species have isotropic Maxwellian speed distributions, and energy is 

63 transferred between species via Coulomb collisions. The simulation spans a rectangular 

64 domain extending 5 Re from Enceladus in all directions except south, where the simulation 

65 extends to 15 Re (Re = 252 km is the radius of Enceladus). 
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A second population of supra-thermal 'hot' electrons is imposed with a fixed density 
(0.3 cm~^) and temperature (160 eV). Hot electrons near Enceladus have been reported 
by Tokar et al. [2009] and have been observed throughout the torus by CAPS and the 
Radio and Plasma Wave Science Instrument [Moncuquet et al, 2005; Young et al, 2005]. 
We showed in Fleshman et al [2009] that a small amount of hot electrons is necessary to 
obtain the ambient ionization. Here we investigate the importance of hot electrons near 
Enceladus itself. 

Neutral source: Following SOS, the plume is prescribed as 



no V ^ 



exp 



e V /r - Rp 



Hq ) \ H, 



(1) 



where Hq = 12° and = 948 km (4 x the Hill radius). SOS offset the plume from 
Enceladus' southern pole by 8° and considered more than one source with the form of 
(1). We consider a single source whose origin coincides with Enceladus' south pole. In 
the nominal case, no is set to 2.5 x 10^ cm~^, corresponding to a neutral source rate of 
fa 200 kg s~^ (SOS). SOS found a much stronger source for EO, so we also investigate the 
implications of a source with no = 2.2 x 10^° cm~^, corresponding to a neutral source rate 
of f» 1600 kg s-^ 

Plasma flow field: Because of the low Alfven Mach number at Enceladus (Ma ~ 
0.1, sutler et al [2008]), perturbations travel rapidly along the magnetic field so that the 
source region presents a cylindrical obstacle to the corotating plasma. We adopt the fiow 
field used by Dols et al [2008] to study the plasma interaction with Jupiter's moon lo: 
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^amb 



cos(20) 
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87 where Vamb ~ 0.8 x Vcor {vcor ~ 26 km s~^) is the ambient plasma speed far from Enceladus 

88 [Wilson et al., 2009]. The magnetic field defines the 2;-axis and is measured from the 

89 flow direction. 

90 Along each flow line, parcels of plasma were started 5 Re upstream of Enceladus with 

91 the steady-state composition given in Fleshman et al. [2009] . The plasma was moved in 

92 the direction of the plasma flow to 5 Re downstream, and the chemistry was updated 

93 at associated time steps. The pickup energy — determined by the relative speed between 

94 neutrals and plasma flow — was also updated. Ions far from Enceladus are picked up at 

95 Vamb-i while thosc near Enceladus are picked up more slowly up- and downstream and more 

96 rapidly on the flanks. 

97 To investigate the effect of ion pickup, we slowed the flow near the obstacle. We followed 

98 Dols et al. [2008] by decreasing the component of the plasma velocity in the flow direction: 

99 itx was replaced by ^{p)u^, where 7(1 Re) = 0.5, increasing linearly to 7(2 Re) = 1. We 

100 find however that a stronger slowing factor (7(1 Re) = 0.1) does not qualitatively change 

101 our results. In Section 3, we compare the pickup rates both for when the plasma has 

102 (nominal case) and has not been slowed. 

103 Two effects are due directly to the slowing of the flow. First, the pickup energy is 

104 reduced, affecting the plasma temperature because fresh ions are picked up at the local 

105 plasma speed. Second, impact ionization increases because plasma transits more slowly. 

106 Impact ionization contributes directly to pickup, as well as indirectly, by seeding multiple 

107 charge exchanges. 
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los We neglect gyromotion on the basis of scale. For example, an H3O"'" ion picked up at Vamb 

109 has a gyroradius of only 0.1 Re- More important is that (in the frame of Enceladus) ions 

no oscillate between zero and twice their pickup speed and we ignore the velocity dependence 

111 of charge exchange. Including gyromotion would enhance the effects we report in this 

112 paper. 

113 Hot electrons: We estimate that hot electrons cool rapidly near the dense plume via 

114 impact ionization, and thus imposed a discontinuity at p = 3Re: 

// -3n Jo 1 < p/Re < 3 

116 However, Dols et al. [2008] showed that field-aligned electron beams (perhaps associated 

117 with the lo auroral footprint) are necessary to model the high plasma density in lo's 

118 wake. At Enceladus, hot electrons may also be related to weak UV auroral spots recently 

119 observed by UVIS (W. Pryor, personal comm.). To investigate the implication of hot 

120 electrons at Enceladus, we consider the additional case where is held at 0.3 cm~^ 

121 throughout the simulation domain. The pickup rates for each case are compared in Section 

122 3. 

123 Pickup rate calculation: Fresh ions are added to the magnetosphere by both charge 

124 exchange (pexch) and impact/photoionization (pioz): 

Aoz = ^ n^UjmjKY^ + ^ nfcmfcK^^°* (4) 

3 k 

126 

Pexch = J]r^('^r^f)m,■/^f'=^ (5) 

i 

128 The reaction rates {Fleshman et al. [2009]) are represented by fi;™^, k'^'^'^^ [cm ^ s~^], and 

129 kF^°^ [s~^]; the ion masses by m; and n^^\ rS^^ are the charge-exchanging neutral and ion 

130 densities. Summations are carried out over processes involving the creation of fresh ions. 
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131 We calculated the time-averages of (4) and (5) to find the average pickup rates for 

132 plasma parcels migrating along each flow line and multiplied by the flow line volume 

133 to find the pickup rate for each flow line. The total pickup rate for each process was 

134 determined by summing the contribution from all flow lines throughout the simulation. 

3. Results 

135 We consider the x-y plane 7.5 Re south of the center of Enceladus. For this, the nominal 

136 case, hot electrons exist throughout the domain, the neutral source rate is Ri200 kg s~^, 

137 and the plasma is slowed in the flow direction (Section 2). The pickup energy, heavy-to- 

138 hght ion abundance (W+/H+), HsO+ZW"*", and HaO"*" temperature (TH30+) in this plane 

139 are shown in Figure 1. The obstacle is plotted in black and 18 flow lines are over-plotted 

140 and labeled. Composition and temperatures along this line are plotted in Figure 2, where 

141 the corresponding flow lines are indicated. 

142 The HsO"*" ion is the most abundant ion in the wake — a consequence explained by 

143 the importance of charge exchange and the fact that charge exchanges lead to HaO"*" 

144 in the presence of an abundant water source. All charge exchanges with H2O in our 

145 model ultimately lead to either H3O+ or H2O+ by H2O+ + H2O H2O + H2O+ and 

146 H2O+ -I- H2O — > OH -I- H3O+. The former reaction supports H2O"'" density somewhat, 

147 but the H2O"'" products also feed into the latter, producing HsO"*". The W+/H+ ratio 

148 increases rapidly because protons are removed by H+ -|- H2O ^ H -|- H2O''". The increase 

149 in electron density is due mainly to impact ionization of the plume by hot electrons. 

150 In Figure 2, the electron temperature has been normalized to the ambient electron 

151 temperature (2 eV), and the ion temperatures have been normalized to their respective 
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152 ambient pickup energies (1.5 and 29 eV). Protons and H3O+ bear the signature of the 

153 cooler pickup temperature from where they were created by charge exchange. A factor of 

154 a few decrease in ion temperature through Enceladus' wake has recently been observed 

155 by Tokar et al. [2009]. The electron temperature has also been cooled by a factor of 2. 

156 Pickup rate: We calculated individually (Section 2) the mass of fresh ions picked 

157 up via charge exchange, Mexch, and Mioz— the ratio of which illustrates the importance 

158 of charge exchange over impact ionization. The simulation was run for the eight cases 

159 shown in Table 1. The flow-field, hot-electron, and source-rate treatments are described 

160 in Section 2. We discuss the cases corresponding to a 'weak' 200 kg s~^ source (Cases 

161 la-4a) and a 'strong' 1600 kg s^^ source (Cases lb-4b) separately. 

162 Weak source ( Cases la-4a): When hot electrons exist locally (la/2a), the total pickup is 

163 roughly 0.3 kg s~^ with a 40% increase when the plasma is slowed. Because of the longer 

164 occupation time associated with the slowed fiow, hot electrons increase seed ionization 

165 (Mioz), and in turn increase charge exchange (Mexch)- When hot electrons are removed 

166 (3a/4a), Mioz decreases by a factor of 4, but Mexch drops by f=:i30%, implying that much 

167 of the pickup is occurring outside the cut-off point at 3 Re. Slowing the plasma has less 

168 effect in 3a than in la because longer occupation does not boost seed ionization without 

169 hot electrons. 

170 Strong source (Cases lb-4b): The effect of the hot electrons becomes more apparent 

171 with the strong source because a denser portion of the plume is intersected by the flow 

172 lines. The effect of slowing the plasma is similar to that in the weak-source case. The total 

173 pickup however is increased by a factor of 3 when comparing the cases with hot electrons 
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174 (lb/2b) to cases without (3b/4b). The almost hnear response of the total pickup to Mioz 

175 (compare lb/2b to 3b/4b) suggests that, unlike in the weak source case, most of the 

176 pickup is occurring inside the hot electron cutoff at 3 Re- 

177 K07 and SOS were in rough agreement on the total pickup rate. For El and E2, they 

178 found ~ 0.2-0.6 kg s~^, and for EO they found ~ 3 kg s~^. Because our model relies on 

179 a physical chemistry calculation alone, it is remarkable to have obtained the same pickup 

180 rates using neutral plume distributions similar to those in SOS (EO: strong source. El and 

181 E2: weak source). 

182 Charge exchange dominates the chemistry by at least a factor of 6 in all cases. Because 

183 of this, the water-group composition ratios (shown only for the nominal case) are quali- 

184 tatively unaffected. In particular, H3O+/W'*" always increases while 0''"/W+, OH+/W+, 

185 and H2O+/W+ always decrease in Enceladus' wake. The dominance of H3O+ elsewhere 

186 has been observed by CAPS during Cassini's orbital insertion period [Tokar et al., 2006; 

187 sutler et al., 200S]. Fleshman et al. [2009] found that a steady-state, water-based Ence- 

188 ladus torus underestimates the H3O+ abundance seen in the CAPS data. Though few new 

189 ions are produced at Enceladus, the process that produces HaO"*" may have an important 

190 effect on the large-scale torus composition. A more complete global model of the torus 

191 should include the effect of dense H2O on Saturn's corotating plasma demonstrated in 

192 this paper. 

4. Discussion and conclusions 

193 To investigate the impact of hot electrons on the chemistry, we have chosen the simplest 

194 flow field possible with minimal perturbation by Enceladus. This flow is roughly consis- 
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195 tent with the compact source derived from magnetometer data (K07) but is much less 

196 perturbed than the flow reported by Tokar et al. [2006]. Similarly, we have started with 

197 a single symmetric plume oriented due-south. In future studies we will explore multiple 

198 jets (SOS), displaced sources (K07), and a minor spherically-symmetric global component 

199 [Burger et al, 2007]. 

200 Our findings are summarized below: 

201 1. Charge exchange dominates the plume-plasma chemistry, confirming previous work 

202 by Burger et al. [2007] and consistent with estimates by Johnson et al. [2006]. 

203 2. Charge exchange leads largely to an HsO^-dominated wake, consistent with INMS 

204 [Cravens et al.., 2009]. Reactions leading to HaO"*" are well known in the comet community 

205 [Aikin, 1974; Haherli et al, 1997]. 

206 3. Comparing our pickup rates to those derived from the Cassini magnetometer (K07 

207 and SOS), INMS, and UVIS [Burger et al., 2007], we find that a persistent source of hot 

208 electrons may exist near Enceladus. If present, beams of hot electrons at Enceladus may 

209 be related to the weak UV auroral spots recently observed by W. Pryor (personal comm.). 

210 Acknowledgments. This work was supported under NESSF grant Planet09F-0036 

211 and CDAP grant 06-CASS06-0062. 
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Case 




kgs ^ 


-^exch 

(10^5 H2O s-1) 


kgs ^ 


(10^5 H2O s-i) 


Mexch/Mioz 


(la) 


Hot electrons + 


slowed flow 


0.25 


(0.84) 


0.038 


(0.13) 


6.7 


(2a) 


Hot electrons + un- slowed flow 


0.22 


(0.74) 


0.034 


(0.11) 


6.6 


(3a) 


No hot electrons + 


slowed flow 


0.18 


(0.61) 


0.0091 


(0.030) 


20 


(4a) 


No hot electrons + 


un-slowed flow 


0.17 


(0.58) 


0.0087 


(0.029) 


20 


(lb) 
(2b) 
(3b) 
(4b) 


Hot electrons + 
Hot electrons + 

No hot electrons + 
No hot electrons + 


slowed flow 
un-slowed flow 

slowed flow 
un-slowed flow 


2.2 
1.9 
0.77 
0.77 


(7.3) 
(6.5) 
(2.6) 
(2.6) 


0.33 
0.30 
0.079 
0.076 


(1.1) 
(0.99) 
(0.26) 
(0.25) 


7.0 
6.6 

9.8 
10 



Table 1. Pickup rates from charge exchange (Mexch) and impact/photoionization {Mioz) for the 



eight cases discussed in the text. Only Mioz increases the plasma density (ne). The cases labeled 
with 'a' correspond to a neutral source rate of 200 kg s~^; those labeled with 'b' correspond to 
a neutral source rate of 1600 kg s~^. Case (la), in bold, is the nominal case from which Figures 
1 and 2 have been generated. 
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Figure 1. Upper-left: the x component of the assumed local plasma flow speed (Section 2) 
normalized to the ambient flow speed (80% of rigid corotation). Upper-right: water-group to 
proton abundance ratio. Lower-left: HaO'^'/W"'" abundance ratio. Lower-right: HsO^ tempera- 
ture. The plane represented here is 7.5 Re south of Enceladus. Model output along the dashed 
line is given in Figure 2. 
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Figure 2. Top: Abundances and electron density from the simulation along the dashed line in 
Figure 1. Bottom: Electron, proton and H3O+ temperatures along the same cut. The electron 
temperature is normalized to its ambient temperature (2 eV), and the ion temperatures are 
normalized to their ambient pickup energies — 1.5 and 29 eV for H+ and H3O"'" respectively. 
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